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A SPECTRAL DECOMPOSITION FOR THE BOLTHAUSEN-SZNITMAN 
COALESCENT AND THE KINGMAN COALESCENT 

JONAS KUKLA AND HELMUT H. FITTERS 


Abstract. We consider both the Bolthausen-Sznitman and the Kingman coalescent restricted 
to the partitions of {1,... ,ri.}. Spectral decompositions of the corresponding generators are 
derived. As an application we obtain a formula for the Green’s functions and a short deriva¬ 
tion of the well-known formula for the transition probabilities of the Bolthausen-Sznitman 
coalescent. 


1. Introduction 

An exchangeable coalescent process is a discrete or continuous-time Markov chain that encodes 
the dynamics of particles grouped into so-called blocks. As time passes, only mergers of some 
blocks may occur, and the rate at which a merger happens only depends on the current number 
of blocks, but not, for instance, on their sizes or the specific particles they contain. The theory of 
exchangeable coalescent processes has its origins in the study of genealogies in population genetics, 
culminating in the seminal work of Kingman [S] . In the context of population genetics Sagitov 
m and later Sagitov and Mohle m derived exchangeable coalescents as limiting genealogies 
of so-called Cannings models, which are discrete-time models of neutral haploid populations 
with exchangeable family sizes. This derivation is also implicit in the work of Donnelly and 
Kurtz [7]. Among exchangeable coalescent processes the so-called A-coalescents have received 
increasing attention in recent years. The latter were introduced independently by Donnelly and 
Kurtz [7], Pitman El and Sagitov m- A A-coalescent > 0} is a time-homogeneous 

exchangeable coalescent process in continuous time with state space Vn, the set of partitions 
of the non-negative integers N := {1,2,...}, that only allows for one merger of blocks at any 
jump. It can be characterized via its restrictions {7T"(t),t > 0} to [n] := (1,... ,n} as follows. 
If at any given time contains b > 2 blocks, then any 2 < fc < 6 of these blocks merge 

at rate Xb^k ■= /g A{dx), where A denotes a finite measure on the unit interval. 

This measure A together with the initial state IT(0) uniquely determines 7T, hence the name 
A-coalescent. 

In this note we consider both the Kingman coalescent > 0} and the 

Bolthausen-Sznitman coalescent — {n^^{t),t > 0}, which are both A-coalescents. For 

convenience we drop the superscripts K and BS when there is no risk of ambiguity. From its in¬ 
troduction Kingman’s coalescent has been used in population genetics as a model approximating 
the genealogy of a sample drawn from a large neutral population of haploid individuals, i.e. in this 
context the particles are interpreted as individuals in the population. The Bolthausen-Sznitman 
coalescent was discovered by Bolthausen and Sznitman [4] in the context of the Sherrington- 
Kirkpatrick model for spin glasses in statistical physics. Goldschmidt and Martin [S] gave a 
construction of via a cutting of a random recursive tree. Bertoin and Le Gall [2] derived 
the Bolthausen-Sznitman coalescent as the genealogy of a continuous-state branching process. 
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An increasing number of recent works suggest that the class of beta coalescents, which con¬ 
tains the Bolthausen-Sznitman coalescent, should provide better null models than Kingman’s 
coalescent for the genealogy of highly fecund populations, see the introduction in [3]. Neher 
and Hallatschek m argue that the Bolthausen-Sznitman coalescent occurs naturally as the ge¬ 
nealogy in models of rapidly adapting asexual populations. Motivated by models of populations 
undergoing natural selection 0, [S], Berestycki, Berestycki and Schweinsberg [T] showed that the 
genealogy of a population governed by Branching Brownian Motion with absorption is again the 
Bolthausen-Sznitman coalescent. 

The article is organized as follows. Our main results are spectral decompositions of the 
generator of the Bolthausen-Sznitman n-coalescent, Theorem [U respectively of the generator 
of Kingman’s n-coalescent, Theorem [2J As Corollaries we obtain for the Bolthausen-Sznitman 
coalescent a derivation of the formula for the transition probabilites that goes back to [3], and a 
formula for its Green’s matrix. As a further application we obtain a spectral decomposition of 
the generator of the block counting process of the Bolthausen-Sznitman, respectively Kingman’s 
coalescent. 


2. Results 


Let us introduce some notation. A partition of a set A is a set, tt say, of nonempty pairwise 
disjoint subsets of A whose union is A. The members of tt are called the blocks of tt. Let #A 
denote the cardinality of A and let Va denote the set of partitions of A. 


2.1. Bolthausen-Sznitman n-coalescent. The Bolthausen-Sznitman n-coalescent = 

^]jn,BS > 0} is obtained by choosing A to be the uniform measure on [0,1], that is the 
corresponding Q-matrix Q := is given by 


( 1 ) 


{ (#p— l)!(#7r—#/9— 1)! 
{#7r-l)! 

— (#7r — 1) 

0 


if TT -< p, 
if TT = p, 

otherwise, 


where tt ^ p if and only if p is obtained by exactly one merger of blocks of tt. In this section we 
only consider the Bolthausen-Sznitman n-coalescent and therefore write TT instead of 

For any two sets A and B C A and a partition tt S Va we call tt\b '■= {C niT: C G tTjCCiB ^ 
0} S Vb the restriction of tt to TT. 


Theorem 1 (Spectral decomposition of the Bolthausen-Sznitman coalescent). Let L = (/,rp) 7 r,pe'P[„] 
and R = (r^p)^_pg-pj^j be matrices defined by 


( 2 ) 


*'7rp 



#7r-#o (*P-1)! 

(#7r-l)! 


ifT^ <P, 
otherwise, 


and 


(3) r^p 




*/ TT < P, 
otherwise, 


where tt < p if and only if each block of tt is contained in a block of p. Then a spectral decompo¬ 
sition of Q is given by Q = RDL, where D = (<Tirp) 7 r,pe'P[„] is defined by (i,rir = —{#tt — 1) and 
d,rp = 0 if TT ^ p. In particular, I^j^tt = = 1 for any tt G V[n]- 


The right eigenvectors r^p may be interpreted as the probability that a random recursive 
tree on the label set tt can be cut down to a tree on the label set p, see the paragraph "A 
connection with random recursive trees" in Section lO below. As an application of this spectral 
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decomposition we derive the well-known formula for the transition probabilities of U given by 
Bolthausen and Sznitman in |3], Proposition 1.4. 

Corollary 1 (Transition probabilities of the Bolthausen-Sznitman coalescent). For any two 
partitions tt, p G V[n] and any time t > 0 the transition probabilities Pjrpit) ■= P{77(t) = p|i7(0) = 
tt} of the Bolthausen-Sznitman n-coalescent are given by 

'' B&p 

where for a; G M, fc G N we denote by '■= x{x -I- 1) ■ • • (t -|- fc — 1) the ascending factorial power 
with the convention x^ := 1. 


Thanks to the spectral decomposition, Theorem [1] we obtain a formula for the Green’s matrix 

= (57rp)7r,pG'P[„] of the Bolthausen-Sznitman n-coalescent defined by Pttp := Jq P-Kp{t)dt. Recall 

that g^p = E[/“ 1 {n(t)=p}dt\II{0) = tt] is the expected total time that 77 spends in p starting 
from TT. We denote by [*] the Stirling permutation numbers or unsigned Stirling numbers of the 
first kind, which count the number of permutations of a set of i elements with j cycles. 


Corollary 2 (Green’s matrix of the Bolthausen-Sznitman coalescent). The Green’s matrix G = 
{gnp) 7 v,peV[„] is given by 


(4) 


g-irp — 


( ^ \ *o 

1 1 (#7r —1)! 


(_l)|fe| 

|fc|-l 


n 


Sep 




00 


0 


i/TT < p {[n]}, 
*/7r < p = {[n]}, 
otherwise, 


where \k\ := Y^Bep^B- 

Remark 1. Notice that since the return probability to any state tt G 'P[n]: t {[s]}, equals 0 
for any R-coalescent 77 there is a close connection between the Green’s matrix G of 77 and its 
hitting probabilities defined by 

7(71, p) := P{77 hits p when started from tt}, 

namely via g^p = h{TT,p)/qp = 7(7r,p)/(l - #p), cf. |I3], p. 146, where qp := 
the total rate in p. 

For a set A and j G N let Vaj denote the set of partitions of A into j blocks. Moreover, {*} 
denotes the Stirling partition numbers or Stirling numbers of the second kind, which count the 
number of partitions into j blocks of a set of i elements. 


Remark 2. From Corollary [5] it follows by a technical but straightforward computation that for 
any tt G and j < i 




pe-PH,: 


(i — 1)! 7 — 1 

k=j 


rJ U 


in agreement with the entries of the Green’s matrix of the block counting process of the Bolthausen- 
Sznitman coalescent provided in the proof of Corollary 1.5 in m- 


As a further application of the spectral decomposition of Q, Theorem [U we derive a spectral 
decomposition of the generator Q' = {qij)i,je[n] of the block counting process {N{t), t > 0} of the 
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Bolthausen-Sznitman n-coalescent defined by N{t) := It is well-known that the matrix 

Q' is given by 

f (i-j)(*-i+i) if * > J) 

<lij = \ if i = j, 

[ 0 otherwise. 

In m this spectral decomposition of Q' was derived by means of generating functions without 
recourse to the partition-valued process U. 


Corollary 3. Let L' = {l'ij)ij^[n], R' = and D' = (d-j)ijeH rnatrices given by 


(5) 



1—1 

1 

i 

bl’ 

1—1 

1 

.11. 

-! 


(*-!)! Ij 

Then a spectral decomposition of Q' is given by Q' = R'D'L'. 


— (f 


2.2. Kingman’s n-coalescent. Kingman ’s n-coalescent 11'^’^ = > 0} is obtained 

by choosing A to be i5o, the Dirac measure in 0, that is the corresponding Q-matrix Q := = 

(g7rp)7r,pej’[„] is given by 


i l if TT < p, 

-{*y) if7r = p, 

0 otherwise, 

where tt < p if and only if tt < p and #tt — #p = 1. In other words, the jump chain of Kingman’s 
n-coalescent is the directed simple random walk on the partition lattice 7^[„], where at each step 
the chain jumps into a coarser partition. From now on we only consider Kingman’s n-coalescent 
and therefore write 11 instead of 77"’^. 


Theorem 2 (Spectral decomposition of Kingman’s coalescent). Define the matrices L = (^^p) 7 r,pePj„] 
and R = (r,rp)7r,pen[„] by 


(7) 



tiT-ttp i*Tr+tp-2)\ 
(2#7r-2)! 


IIbgp 


ifT^^P, 

otherwise, 


and 


(8) r^p 


(2<!p-l) 


(#7r+#p 

0 






ifTr<P, 

otherwise. 


Then a spectral decomposition of Q is given by Q = RDL, where D = (d^pj.n-^pg-pj^j is defined by 
c^TT-n- := ~{* 2 ) and d^rp = 0 if tt y p. In particular, = 1 for any tt € 'P[n]- 

Remark 3. Although the spectral decomposition of Q, Theorem [5J directly yields a spectral 
decomposition of the matrix of transition probabilites p,rpit) ■= P{77(t) = p|7T(0) = tt}, analogu- 
ously to the one given in for the Bolthausen-Sznitman n-coalescent, this expression is not 
particularly handy. In fact, a better approach to calculating pTrp(t) is the one given by Kingman 
in [9] Equation (2.5). 


As in the case of the Bolthausen-Sznitman coalescent, consider the block counting process 
{N{t),t > 0} of Kingman’s n-coalescent defined by N{t) := #n{t) and let Q' = {qij)ij^[n\ 
denote the generator of N{t). From Theorem we obtain a spectral decomposition of Q'. 
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Corollary 4. Let L' = R' = be matrices defined hy 


(9) 


7 / . ( ^)' 

% 

.J . (2J-1)! 

i 

^ ’ (2z - 2)! 

j. 

’ (* + J-l)! 

j_ 


where [*J := denotes the unsigned Lah numbers which count the number of partitions into 

j blocks of a set of i elements, where the elements in each block are ordered. Then Q' = R'D'L' 
is a spectral decomposition of Q', where D' = is defined by dh = —{^) ^ij = 0 */ 

i 7^ j- 


3. Proofs 

In order to prove our results, we need some notions and facts from the theory of lattices, that 
we collect from HZ]. Recall that a partially ordered set (poset for short) (P, <) is a set P together 
with a binary relation < satisfying: 

(1) For all p G P, p < p (reflexivity). 

( 2 ) If p < g and 9 < p, then p = q (antisymmetry). 

(3) If p < g and q <r, then p < r (transitivity). 

Recall that two posets P, Q are called isomorphic, in which case we write P = Q, if there exists 
an order-preserving bijection (f: P ^ Q whose inverse is order-preserving. The cartesian product 
P X Q oi two posets P, Q is defined on the set {{p,q). p € P,q & Q} by letting (p, q) < (p', g') 

in P X Q if and only if p < p' in P and g < g' in Q. For p, g G P an upper bound of p and 

g is an element r G P such that p,q < r. A least upper bound of p and g is an upper bound 
s € P oi p and g such that for any upper bound r of p and g one has s < r. Clearly, if a least 
upper bound of two elements p and g exists, it is unique. A greatest lower bound is defined in 
complete analogy. A lattice P is a poset with the property that any two of its elements have 
a least upper bound and a greatest lower bound. It is well-known, that P[n] together with the 
relation < as defined in Theorem [T] is a lattice, the so-called partition lattice. For tt, p G P[„] 
with TT < p we call the set [tt, p] := {a G P[n] : tt < ct < p} an interval. We will make repeated 
use of the isomorphism 

( 10 ) [7r,p]^ X IP.lB. 

Bep 

cf. Example 3.10.4 in (TT]. For more information on posets in general and the partition lattice 
in particular the reader is referred to HZ). Evidently, we have V[n] = {W}]) where for any 

set A we let A a '■= {{«}: a G A} be the partition of A into singletons. 

Using the notation we just introduced, the n-A-coalescent TT" is a Markov chain with state 
space T’fn]; initial state Z\[„] and Q-matrix Q" = (g"p) 7 r,pg'P[„] given by 

{ '^#'7r,#7r—#p-t-l if TT ^ p, 

-A#^ if TT = p, 

0 otherwise, 

where Af, := Y.l =2 (D^b.k, b>2, is the infinitesimal rate. There are several extensions of < to 
a linear order on R[„]. Let us fix such an extension <ex and notice that the following quantities 
of Q" that we are interested in do not depend on the specific extension chosen. The linear order 
<ex induces a natural bijection tp from P[n] to [Bn], where Bn denotes the nth Bell number, 
which is the number #R[„] of partitions of [n], defined inductively by letting '!/i(Z\[„]) = I and 
1 / 1 ( 71 ) < '(/’(p) iff TT <ex p. Then Q" can be seen as an upper right triangular matrix with entries 
ordered according to <ex, if we define row/column tt to be lower than row/column p iff tt <ex p. 
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The determinant of Q" is therefore given by the product of its diagonal entries. Hence, the 
characteristic polynomial of Q" is given by 


n 


XQn(a:) := det(Q” - a;I„) = (-1)^” (^#^ + 2^) = (-1)^" 


where is the identity matrix on 'P[n]- Hence, for each is [n], —Xi is an eigenvalue of Q" with 
algebraic multiplicity #V[n],i = {"}• Since Q" is a Q-matrix, ( 1 ,..., 1 )^ S is an eigenvector 
corresponding to the eigenvalue 0, and direct inspection yields that (1,0,S is also 
an eigenvector corresponding to the eigenvalue —A„. From now on we fix an n S N, n > 2, and 
drop this index in the notation, if there is no risk of confusion. 


3.1. Bolthausen-Sznitman n-coalescent. In order to prepare the proof of the spectral de¬ 
composition of Q = TheoremlU we calculate the left and right eigenvectors of Q. In the 

sequel we give two proofs for the right eigenvectors of Q that are of rather different flavours. The 
first proof is completeley self-contained and only makes use of the partition lattice ■ Together 
with the proof of Lemma [S] it might serve as a starting point to find a spectral decomposition 
for more general coalescents, e.g. beta coalescents. There is a probabilistic interpretation of the 
right eigenvector of Q in terms of random recursive trees which then motivates our second proof 
that heavily draws on random recursive trees and their connection to the Bolthausen-Sznitman 
coalescent as explored in Goldschmidt and Martin [S]. 

Lemma 1. For p £ 'P[n\ the vector (r.n.p),n.g-pj^j defined by 



( 12 ) 


is a right eigenvector of Q with corresponding eigenvalue 1 — #p. 

Proof, (first proof of Lemma [1} In order to carry out the following calculations, for any two 
partitions tt, p £ '^[n] "we need to parameterize the set {a: n ^ a < p}. To construct an arbitrary 
partition a such that tt ^ cr, we could choose a subset C C tt of at least two blocks of tt and 
merge them in order to obtain a. In this case #a = #tt — #C -1-1. If, additionally, we require 
a < p, we certainly cannot choose any collection C of blocks in tt. Instead, all blocks chosen have 
to be in ■k\b for some block i? £ p, in which case #(t|s = #tt\b — #C + 1. To summarize, we have 
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Using this parametrization, we obtain 

= X! ^TvaTcrp “ (#71 “ ^)r^p 

(7 : 7r^cr<p 

(#7r-#C')!(#C-2)! (#p-l)! 


(#77 — 1)! (#77 — #C)! 


= E E 

-BGpCCti-Ib 
#C>2 

X (#77|s - #(7)! (#77|£) - 1)! - (#77 - l)r^p 

Dep\{s} 

- E i: n<*'b - c - c - d. 


BSp c=2 ^ 

#7r|B 


' (#77|b - 1)! 


TTp 


D€p 


= I E E 2)!7S^1) I’■•p 

. BGp c=2 ^ ^ 


' {#7t\b - 1)! 


#7r|B 


= E 


1 


^ c(c — 1) 

BGp c=2 ^ 7 


- (#77 - 1) 


= ^ (# 77 |s - 1 ) - (#77 - 1 ) I r^p = (1 - #p)r^p, 

\Bep 

and the claim follows. 


□ 


A connection with random recursive trees. Let us now recall the notion of a random 
recursive tree in order to prepare our second proof of Lemma [1] where we closely follow Gold¬ 
schmidt and Martin [H]. We call a tree on n nodes labelled by 1,2,... ,n an increasing tree if 
the root has label 1 and the labels in any path from the root to another node are increasing. 
Figure [T] shows all 6 increasing trees on {{1}, {2}, {3}, {4}}. If we have an increasing tree on n — 1 


{4} 



Figure 1. All 6 increasing trees on {{!}, {2}, {3}, {4}}. Only the third and the 
sixth tree contain {{1, 2, 3}, {4}}. 

nodes, we obtain an increasing tree on n nodes by adding a node labelled n and attaching it by 
an edge to one of the nodes in the given tree. Starting from the tree that only consists of the 
root node 1 this gives an explicit construction of all increasing trees on n nodes. Consequently, 
there are (n — 1)! increasing trees on n nodes. A random recursive tree is a tree chosen uniformly 
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at random from all increasing trees on n nodes. An explicit construction of a random recursive 
tree on n nodes is the following. Start with the root node labelled 1. If the tree has k nodes, 
choose one of these nodes uniformly at random and attach to it node fc + 1 by an edge. Stop 
after attaching node n. 

For any partition tt G 'P[n\ sn increasing tree on tt is a tree with #tt nodes that are labelled 
by the blocks in tt such that the labels in any path from the root to another node are increasing 
with respect to their least element. We denote a random recursive tree on tt by Ttt. 

Crucial to the construction of the Bolthausen-Sznitman coalescent via random recursive trees 
is the following cutting procedure. When given a tree 7^ on tt, a cut is performed by picking an 
edge, removing the subtree above this edge (here we picture trees as they grow in nature: from 
the root at the bottom to the leaves at the top) and adding the labels of this subtree to the labels 
of the node below the edge. For a simple tree the cutting procedure is depicted in Figured! We 


{4} 



Figure 2. An increasing tree on {{!}, {2}, {3}, {4}} cut down successively to 
a tree on one node {1, 2, 3,4}. 

denote by cT the tree obtained by cutting T at an edge chosen uniformly at random. A striking 
result then is Proposition 2.1 in Goldschmidt and Martin [H], which states that after cutting a 
random recursive tree on tt at an edge chosen uniformly at random, the new tree is a 
random recursive tree on the new label set (which is again a partition of [n]). 

For TT G Vin] define a time-homogeneous continuous-time Markov chain 77^ := {7?.7r(t),t > 0} 
with values in the set of random recursive trees as follows. The initial state 7?.7r(0) is the random 
recursive tree T,r on tt. The process 7^^ evolves according to the following dynamics. Suppose 
TZtt is in state T- If T has only one vertex, do nothing. Otherwise, attach to each edge in T 
an exponential 1 clock, all clocks being independent. When the first clock rings, cut T at the 
associated edge to obtain the next state of TZj^. For any tree T on the label set tt let piT) = 
TT. Proposition 2.2 of Goldschmidt and Martin [S] establishes that the Bolthausen-Sznitman n- 
coalescent U is equal in distribution to 

{p(nA„{t)),t > 0 }. 

Fix two partitions tt, p G 'P[n] such that tt < p. We say that an increasing tree T on tt contains 
p iff one can obtain an increasing tree on p by successively cutting T'. If T' contains p we write 
p \ZT. Let {RTr{k),k > 0} denote the jump chain of 77^. From the definition of TZ it follows 
that a process that is equal in distribution to the jump chain R can be constructed recursively 
by letting 

R-nifi) '■= ‘^TT, R-jvik + 1 ) := cRTr{k), k > 0. 

A simple question then is: what is the probability P{p C that a random recursive tree on 
TT contains pi Notice first, that an increasing tree T on tt contains p if for each block B G p 
we can find a node v = v{B) in B such that the labels in the subtree above v coincide with the 
elements in ttI^. In other words, we can construct all increasing trees on tt that contain p by first 
constructing an increasing tree T on p. Then each node f G T is labelled by some block B G p. 
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Now if #7r|B > 1, build an increasing tree on ttI^ and replace v by this tree. This procedure is 
done for all nodes v G T to obtain an increasing tree on tt that contains p. Therefore, the number 
of increasing trees on tt containing p is 

#{T: T increasing tree on 7r,p C T} = (#p — 1)! i#Tr\B — !)!• 

B£p 


On the other hand, the total number of increasing trees on tt is (#7r — 1)!. Since a random 
recursive tree on tt is a tree chosen uniformly at random from all increasing trees on tt, we obtain 


(13) F{p □ TJ 


(*p-i)!nBgp(*^iB-i)! 

(#7r —1)! 


if TT < p, 

otherwise, 


which is just r^^p. We can now turn to our second proof of Lemma[T]in terms of random recursive 
trees. 


Proof, (second proof of Lemma 1) We rewrite the statement as 
(14) ^ q^^arap = (#7r - #p)r.^p. 

<7 : TT^cr 

Let J = {J(fc),fc > 0} denote the jump chain of 77 = Since P{J(1) = (t|J(0) = tt} = 

is the probability that J jumps from tt to a, after dividing (|14l) by = #7r — 1 we obtain 
for the left hand side 

E —E P{ J(1) = al J(0) = 7r}P{p C T4 

a: TT^a a: tt^ct 

= E P{p(cT,) = a}P{p c ^4 

( 7 : 7r-<cr 

= P{pETp(,5j} 

= P{p c cT^} 

= P{p c cT^Ip □ T,r}P{p c 'Itt}, 


where we used that for any random recursive tree P one has —d P- Conditional on p C Ttt, 

there are precisely #7r — #p among the #tt — 1 edges in that we may cut in order to obtain a 
tree that contains p. Since, by definition, cT,r is obtained by cutting at an edge chosen 
uniformly at random, we have 

P{p c cT^Ip c 

#71 — i 

The claim follows. □ 


Lemma 2. For a; G R \ {0} we have 

(15) E n (-4^- 

crGP[„] ^ '' Bep 

Proof. Notice that 

{7GV[n] Cr£[TT,p] B^(7 

^ ’’ a^[Tr,p\B^a 


(16) 
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From the interpretation of [”] as counting the number of permutations of [n] with i cycles it is 
clear that ["] = ^ ~ 1)^ equation (1.15) in [H]. Using the isomorphism 

[tt, p] = X^gp 'Ptt\b I ■''^e calculate 


(7^[k,p] B^cr 


(17) 

(18) 


= n E n-"(*e-i)! 

BepreV^lg CGt 

( #7r|B 

E(-)‘ E n(‘e-i)! 

k = l Te'P,r|B,fc 

#7r|B 

= nE(-)‘ 

B^pk—1 

= 

Bep 


#tt\b 

k 


where in the last step we used x” = Y^=i ef. equation (1.16) in [T?| . 

Lemma 3. The matrix L = {lTTp)-n,p^'P[n,] defined by 

(19) i. 

^ ^0 otherwise, 

is the inverse matrix of R, i.e. X]crG'P[ ] 

Proof. Choosing a: = 1 in Lemma we have that X]itgP[ ] ^tto'^ctp = ^^^p- 
Proof, (of Theorem [1]) The claim follows by Lemmata [1] and |3l 


□ 


□ 

□ 


Proof, (of Corollary [T]) Since 77 is a continuous-time Markov chain with finite state space, we 
have for P{t) = (p7rp(7))-7r,pe7:’[„] the identity P{t) = exp{tRDL) = Rexp{tD)L. In the last step 
we made use of the spectral decomposition, Theorem 1. In particular, this yields 


(20) p^p{f)= ^ r^tre ^"'lap= ^ '77<T 

<y^'P[n] G^[k,p] 

Letting x = e“* in Lemma proves the claim. 
Proof, (of Corollary [21) From ([T71) we have 


t(#cr-l)7 

' Trrr^ ^ap' 


□ 


#7I-|b 




X * ^ P[ -a:(#7r|B - 1)! = x ^ H LI 

cre[7r,p] BGct B^p fc—1 


#7r|B 

k 


= E (-1) 

(feB)Bep6N»P 




n 

BGp 


#711 s' 
kB 


where |fc| := X^sep ^b- Letting x = e * and integrating out with respect to t we see for p ^ {[n]} 


E Y[-e \#T^\B-iy-dt= 


(_l)|fc| 


n 


crG[7r,p] B^a 


\k\ -1 

(feB)Bep6N»P ' ' Sep 


#7r|B 

kB 
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where for p = {[n]} 

pOO 

o-efTi-.p] Bea k=i 

The claim follows from (USD and the definition of g^^p- D 

Proof, (of Corollary [SD Evidently, the rate at which a jump from i to j blocks occurs equals the 
sum of all rates at which a jump from a partition tt G 'P[n],i to ^ partition p € 'P[n\,j occurs, 
i.e. q'ij = 971-p, where tt e T’[„]y is fixed arbitrarily. The quantity 

not depend on the choice of tt G the following calculation shows. By the spectral 

decomposition of Q, Theorem [TJ we obtain 




#7r 

k 


= E E 

P^'P[-n],j (7G[7r,p] 




^ (#0' - 1)! 

(#7r — 1)! 

(7: 7r<crp: (7<p ' B^cr 

#(T>j #p^j 




(j-l)! 


#<T-j 


(#a — 1 )! 


(J-I)! , 

(*-!)! 


f—11* 

- 1 )^ ^ 1 ) 1(1 


<7: 7r<(7 
#tT>j 


(*-!)! 
hence, = R'D'L'. 








(i-fc) 


#cr 


(-1)^ 
(fc-l)!b- 


(#CT - 1)! ( j 


^ ^ ^ik^kk^kj 5 
k—i 


□ 


3.2. Kingman’s n-coalescent. For any two partitions tt, p G V[n] such that tt < p one may ask 
in how many different ways the jump chain of Kingman’s coalescent may reach p when started 
in TT. Since at each step only one merger of a pair of blocks occurs, there are #7r — #p + 1 steps to 
be taken, and so the set of different ways is C{tt, p) := {(tti, ..., tt^) : tt = tti < • • • < = p,m = 

#7r — #p + 1}, where we defined < in the paragraph preceding Theorem [H We call each element 
in C{tt,p) a maximal chain in [tt, p] and denote by 771 ( 71 , p) := #C(7r,p) the number of maximal 
chains in [7r,p]. Before we turn to the proof of the spectral decomposition of Q, Theorem [5J we 
count the number of maximal chains m(7r, p) in [tt, p] in the next Lemma. 

Lemma 4 (Number of maximal chains). For 7r,p G P[n] with tt < p we have that 

(21) 777 ( 71 ,p) = {tfTT — tfp)l ^ttIb!. 

BGp 

Proof. Notice that any maximal chain (tti, ..., 7r„) in [/![«], {[ 77 ]}] can be constructed as follows. 
Let TTi := Z\[n], and if tt^ with z < 77 is constructed, tt^+i is obtained by merging two blocks in tt^, 
which can be done in (* 2 ’) ways. When 7r„ = {[ri]} is reached, the construction is finished. This 
construction shows that there are 777(Z\[„], {[ 77 ]}) = ( 2 ) ("b) ’'' (D “ 2^“"77!(77 — 1)! maximal 
chains in {[ 77 ]}]. Hence (1^ holds in the case (7r,p) = (Z\[„], {[ 77 ]}). 

For the general case, recall the isomorphism [tt, p] = X^gp^Tri^. As a consequence, any max¬ 
imal chain in [7r,p] can be built by choosing a maximal chain in each factor T^ttIb then 
’intertwining’ these chains, i.e. ordering their elements (excluding the first element — the par¬ 
tition into singletons — in each chain) in any order subject to the restriction that the order of 
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elements of the same chain is preserved. Consequently, we have 

m{Tr,p) = (#7r-#p)! [(#7r|B - 1)!]"^ 

BSp B£p 

= (#7r — #p)\ = 2*'’“*^(#7r — #p)\ #tt\b- 

BSp BGp 


□ 


Lemma 5. For any p G 'P[n\ the vector (r 7 rp) 7 re'P[„] defined by 


( 22 ) r^p 

(23) 


(#7r-#p)!(#7r+#p-l)! "H'' ) P! 

0 


otherwise, 




otherwise, 


is a right eigenvector of Q with corresponding eigenvalue —{* 2 )- 


Proof. Fix TT, p G '^[n] • If tj < p, we have 

^ ( Q-TTcrTcrp 

o-e'Pr„i 


= E 


' (7p 


#7r 

2 


•Kp 


-*P{2#p-l)\ 


(#7r - 1 - #p)!(#7r + #p - 2)! 

(#7r — #p)(#7r + #p — 1) f#TT 

2 


E ^A^p)- 


#7r 


' TT/O 


where we used ir<o- ^( 0 ") p) = rn{TT, p). If tt = p, we have 


^ ^ QTZa’Pap - ^TT-TT - 


Trp, 


since 771 ( 71 , 71 ) = 1, hence r^ 7 r = 1- Finally, if 71 < p does not hold, thus r^p = 0, we cannot have 
T < (T < p for any a G P[„] and therefore Q^rarap = 0. This shows (1^ . Now (1^ follows 

from Lemma 0] on the number of maximal chains. □ 


Evidently, the Bn eigenvectors of Q defined by (12^ are linearly independent. We are now 
interested in the inverse matrix of i? = (r7rp)7rpg77[^j, that is the matrix L = {Ittp)tt,p^V[„,] such 
that ^Trp — ^^^Ap fol" ^11 P £ 


Lemma 6. For any ti G P[n] the vector (l,n.p)pg-pj^j given by 

[ t) (2#7r-2)!(#7r-#p)! , p) */71 S P, 


(24) 

(25) 


t-Kp - 


0 otherwise, 

*/7i<p. 


otherwise, 


is a left eigenvector of Q with corresponding eigenvalue —{* 2 )- 
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Proof. Use a<p O') = p) to obtain 

^ ) ^TTcrQcrp 
-TGP[„, 

- V / - I 

— 2_^ 12/^^ 

a: <T<p ^ ^ 


(2#7r — 2)!(#7r — #p — 1)! 

(#7r + #p — l)(#7r — #p) ^#p 

2 


O': cr</3 


^TTp — 


f-TTp 


•7rp5 


thus (I24II holds. Equation (12511 follows from the Lemma on the number of maximal chains, Lemma 

m □ 


Proof, (of Theorem [5]) The inverse matrix of i?, let us call it 17 = , is uniquely 

determined, and is a matrix of left eigenvectors of Q, i.e. UQ = DU. Moreover, for any tt S V\n\ 
we have by assumption 7 r<(T<p ^’ro-Co-p = iJu-tt = 1- This uniquely determines a matrix 

of left eigenvectors of Q, since for any tt, p G P[n] with tt < p we have 

II-Trp — ^ ) Il7r(T^crp — l{7r<p} T ^ ) UTrcrQ^crp: 

o-gf-n-jp] O’: -ir<a<p 

and Uttp = 0 for tt ^ p. Since Qi? = RD by Lemma ([ 6 ]), and evidently = 1 for any tt € P[n]^ 
we have U = L and the claim follows. □ 



Remark 4. Instead of calculating the hitting probabilites h{TT,p) of Kingman’s n-coalescent via 
the spectral decomposition. Theorem [5J we use the observation from section that the jump 
chain of TT can be interpreted as the directed simple random walk on 'P[n\ ■ This implies that the 
jump chain of TT (when started from tt) traces out any maximal chain in [tt, {[n]}] with equal 
probability 771 ( 71 , {[n]})”^. Clearly, the total number of maximal chains in [tt, {[n]}] that contain 
p is m(TT, p)m{p, {[n]}) , and thus 


h{TT,p) 


m{TT,p)m{p,{[n]}) 

w(7r,{[n]}) 


#7r — 1 

#p — 1 



#7r|B! = 

B^p 


#TT 

.*P. 


n 

Bep 


where we used the Lemma on the number of maximal chains. Lemma 01 in the second step. 
In the special case tt = Z\[„] this formula was given by Kingman in |^, equation (2.3). 


Proof, (of Corollary 0]) In complete analogy to the argument in Corollary 01 we have q[j = 
X)p 6 B[ ] Tttp independent of the particular partition tt S P[n],ii the following calculation 
shows. Using the spectral decomposition of Q, Theorem 01 we obtain 


E E ' 

PGT’hj o-GK.p] 


= 


7T<T^<T<T'^(7p 


= - E 

(T: 7r<(T 
#cr>j 


(2#cr - 1)! 
(i + #CT — 1 )! 



/#(t\ / ,#g+j (o' + j-2)! 

I 2 / ^ ^ ^ (2#a-2)! 

^ ' P'- <^^P 

#p^j 


#cr|B! 

BGp 


^ ( 2 fc-l)! 

^ (7 + fc - 1)! 


i 

k 


AV (fc + 7 - 2)! k 
{2j^ > {2k-2)\ Lu 


I 

k=j 
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where in the third step we used the identity ^ risecr twice. This identity is 

obvious from the interpretation of [*J as the number of partitions into k ordered blocks of a set 
of i elements, where the elements in each block are ordered. □ 
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